1. Joint Modeling
1.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
1.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.002 1.01
## LDevsum 0.999 1.00
## dh0 1.918 3.57
## dh1 1.668 3.63
## dh2 1.135 1.45
## dl0 1.658 2.87
## dl1 1.334 2.03
## dl2 1.074 1.29
## dl3 1.258 1.85
##
## Multivariate psrf
##
## 1.49
1.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
1.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1231.31073 | 22756.8074 |
| DIC3 | 1164.66422 | 20662.5287 |
| PWAIC | 46.58249 | 994.8117 |
| WAIC | 1194.22015 | 21039.8413 |
2. Separate Modeling
2.1 Trace plots for convergence check
The current MCMC setting is:
- 300,000 iteration;
- 200,000 burn-in;
- 100 thinning.
2.2 Gelman and Rubin’s convergence check
## Potential scale reduction factors:
##
## Point est. Upper C.I.
## HDevsum 1.00 1.01
## LDevsum 1.00 1.02
## dh0 1.23 1.79
## dh1 1.12 1.40
## dh2 1.02 1.09
## dl0 1.11 1.12
## dl1 1.15 1.51
## dl2 1.03 1.11
## dl3 1.19 1.37
##
## Multivariate psrf
##
## 1.21
2.3 ACF Plots
Here we plotted ACF plots for the following variables:
- Total deviance;
- Variables that didn’t pass the convergence check.
2.4 WAIC results
| LevelH | LevelL | |
|---|---|---|
| DIC | 1377.52079 | 23110.679 |
| DIC3 | 1247.01158 | 20729.655 |
| PWAIC | 90.84391 | 1040.447 |
| WAIC | 1316.35844 | 21133.344 |